Interface Tension of Bose-Einstein Condensates 
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Motivated by recent observations of phase-segregated binary Bose-Einstein condensates, we pro- 
pose a method to calculate the excess energy due to the interface tension of a trapped configuration. 
By this method one should be able to numerically reproduce the experimental data by means of a 
simple Thomas-Fermi approximation, combined with interface excess terms and the Laplace equa- 
tion. Using the Gross-Pitaevskii theory, we find expressions for the interface excesses which are 
accurate in a very broad range of the interspecies and intraspecies interaction parameters. We also 
present finite-temperature corrections to the interface tension which, aside from the regime of weak 
segregation, turn out to be small. 

PACS numbers: 03.75.Hh, 67.60.Bc, 67.85.-d, 67.85. Be, 67.85.Fg 
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I. INTRODUCTION 

Soon after the first experimental realization of Bose- 
Einstein condensation (BEC) in diluteeases, binary BEC 
mixtures were realized [H, S, i, i, S i, 0, i, H ES]- 
These experiments precipitated a strong theoretical in- 
terest, the origin of which is the fact that the multi- 
component BEC is not a trivial extension of the single- 
componcnt BEC 



[30 |. Indeed, novel 
and fundamentally different physics arise on both micro- 
scopic and macroscopic levels. Phenomena studied for 
binary mixtures so far include the quantum tunnelling of 
spin domains [jj_p|, spin-relaxation processes [4|], vortex 
configurations 0, formation of Feshbach molecules @ 
and collective oscillations 0]. Experiments with phase- 
segregated BEC mixtures were established already in 
1998 ^Q, lZ(]; this triggered many theorists to focus on 
the observed weakly-segregated phases whereby the im- 
portance of the interface physics was highlighted several 
times 18, 19, 20, 21, 22, 23, 24, 25, 26, 27]. 

Also surfaces of .smj^Ze-component BEC gases have 
gained much attention. Studies are performed on vor- 
tices, surface modes and tunnelling phenomena at the 
border of a trapped single-component BEC [U [H, [H, 
[3^ and on BEC diffraction and van der Waals forces us- 
ing an optical mirror or evanescent-wave prism [ssl . [sgI . 
|37| . Moreover a combination of the interface physics of 
binary BEC gases and the surface physics near hard walls 
led to the theoretical prediction of anomalous wetting 
phase transitions [ssj . 

The clearest observation so far of phase segregation of 
BEC mixtures was reported by Papp et al. (see Ref. [lol ) 
who used a mixture of ^^Rb and ^^Rb particles. By 
changing the particle numbers of both species and the 
intraspecies interactions of the ^'^Rb particles by use of 
a Feshbach resonance, many topologically distinct states 
were encountered. In the phase-segregated regime most 
trap configurations contained BEC droplets. Papp et 
al. stated that a detailed theoretical understanding of 



the observed droplet formation is lacking. As an essen- 
tial first step towards this understanding we present here 
expressions for the interface tension. Moreover, we also 
formulate a method by which all trap configurations can 
be straightforwardly studied by the use of a local den- 
sity approximation or Thomas-Fermi approximation. At 
first sight, numerically reproducing the experimentally 
observed (meta) stable states requires solving the cou- 
pled Gross-Pitaevskii (CP) equations for the full three- 
dimensional system; in such case, a high accuracy is indis- 
pensable to capture the important energy contributions 
near the two-phase boundaries which strongly affect the 
ground state topology. 

In this work, we argue that it is not necessary to solve 
the full GP system in order to recover numerically the 
experimental configurations; instead, one can find the 
ground state, that is, minimize the total energy of the 
trapped cloud, by use of the Thomas-Fermi approxima- 
tion and the expressions for the interface tension which 
we provide here. This allows an omission of the quantum 
pressure term which strongly complicates the numerical 
work. The analytical calculation of the interface tension 
is nontrivial as it involves solving the binary GP equa- 
tions (i.e. two coupled nonlinear second-order differential 
equations). Estimates and analytical expressions were al- 
ready given in Refs. [2l|j2^ and were relevant to the early 
experiments on weakly-segregated mixtures but irrele- 
vant to the interfaces as observed in Ref. [l^ since there 
segregation is not weak. By developing expansions, we 
succeed in obtaining analytical expressions for the inter- 
face tension which are accurate for almost all parameter 
regimes. 

As a main result we find that, even though the interface 
tension is not constant throughout the trap, the total 
excess energy fi^ can be reduced to a simple integral of 
the trapping potential along the interface area A: 



= / dr[M,-C/,(r)]3/2, (1) 

where a^^ , , /i^ and [/^ are the intraspecies scattering 
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length, the particle mass, the chemical potential and the 
applied trapping potential of phase 1. The parameters K 
and ^2 /^^ can be expressed in terms of the atomic masses 
and the interspecies and intraspecies scattering lengths: 



K 



"^-+'"- and4/C.- 



(2) 



For a broad range of values of /^^ and K, the function 
can be written as: 



HU(.,K)^^a + UO-"-^^^ (3) 

/0.055 0.067 



This work is structured as follows. We start off in 



Sect. [IT] by introducing the GP formalism for binary 
BECs. After defining the excess energy in a homoge- 
neous and in a trapped system in Sect. Illli we find an- 
alytical expressions for the interface tension in Sect. IIVI 
We estimate the finite-temperature corrections to the in- 
terface tension in Sect. |Vl In Sect. IVII we discuss the 
experimental relevance of our findings. 



II. BINARY BEC SYSTEM 

Equilibrium states of a mixture of two dilute BEC 
gases with order parameters ^'^ and 4" 2 and chemical 
potentials /x^ and /x^, are well modelled using the grand 
potential: 



%l*.(r)|^ 



/ dr|*,(r)n*2(r)|^ (4) 



where G — 2Tih^a..(m 



^) are the coupling con- wave functions ip-^ and ip^ 



stants and a., are the s-wave scattering lengths (hence- 
forth i = 1, 2). In view of the derivation of the interface 
tension, we continue here by assuming vanishing external 
potentials U-^ — — 0; we reintroduce the external po- 
tentials when discussing the surface excess of a trapped 
system. In the absence of particle flow one can choose the 
order parameters to be real valued such that the equilib- 
rium pressures for the pure states are: 



P 



2G,. 



(5) 



In order to study phase-segregated states of the two co- 
existing condensates, we introduce the parameter 



K I 



(6) 



which quantifies the interspecies couplings relative to the 
average of the intraspecies couplings. It is well-known 
that when K exceeds one, the species become immiscible 
and only pure phases can exist [2^ while phase mixing 
occurs when K < 1. Along the two-phase interface, co- 
existence is ensured by the condition: 



A. 

2G,, 



2G, 



(7) 



Henceforth we assume that K > 1 and that the chemical 
potentials are such that the bulk pressures are = = 
P. For calculational convenience, we further rescale the 
order parameters ^'^ and ^'^ and define the dimensionless 



(8a) 



The quantum nature of the system results in a zero-point 
motion; this determines the typical length scale for den- 
sity modulations at boundaries, impurities, vortices or 
solitons. For the pure phases, the resultant length is the 
healing length which is defined as: 

Without loss of generality, we choose the phases such that 
^2 /^^ < 1 . As coexistence must occur along the interface, 
K and £,2/^1 can be expressed in terms of the atomic 
masses and the scattering lengths as given in Eq. ([2]) . Fi- 
nally, after rescaling the space coordinate to the dimen- 
sionless variable r = r/^^ , the reduced time-independent 
Gross-Pitaevskii (GP) Eqs. are [10,1411: 

VVi =-V'i +^' + i^V'iV'2', (10a) 

= - ^2 + € + ^V-.V-'. (10b) 

Having established the equations of motion, we continue 
now by defining the interface tension. 



III. DEFINITION OF INTERFACE TENSION 
AND INTERFACE EXCESS 

In order to calculate the interface tension, consider two 
BEC components in an infinitely large system with trans- 
lational symmetry in the x-y direction. The presence of 
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phase 1 at z ^ cxD and of phase 2 at z ^ — oo imply the 
following boundary conditions: 

■(/'i(z = — oo) = ijj2{z = oo) = 0, (ll^-) 
■0^(2 = oo) ~ tp2{z ^ — oo) = 1. (lib) 

Under these conditions, the first integral of the reduced 
GP Eqs. pHal) and pOb| is: 



^ V ' 2 

= 1,2 



(12) 



where the dot denotes the derivative with respect to z. 
Since we work at fixed chemical potentials the interface 
tension is the excess grand potential per unit area; this 
excess is uniquely determined by subtraction of the total 
grand potential of a volume V containing a pure phase, 
from the grand potential fl, that is, + PV. By use of 
Eqs. PU)) and the interface tension can be written 
as: 



(13) 



The interface tension in the grand canonical ensemble is 
exactly four times the interface tension in the canonical 
ensemble as introduced by Ao and Chui [l^l [l^- Note 
also that, as opposed to the grand canonical ensemble, 
the definition of the interface tension in the canonical 
ensemble is not unique [s^. In Eq. we define the 

dimensionless function which, as is also the case for 
the normahzed profiles ^p., only depends on ^2/^1 ^ 
since the wave functions are fully determined by the GP 
Eqs. pO|l and the boundary conditions (fTT|) . It is clear 
from Eq. that the excess energy is positive and due 
to a bending of the normalized wave functions V'l and ip^ . 

One can now straightforwardly generalize the defini- 
tion of the interface tension to the total excess energy 
of a trapped system. Therefore we use a local approxi- 
mation for the interface tension: we assume the charac- 
teristic lengths over which the trapping potentials (r) 
(for i = 1, 2) vary to be large as compared to the in- 
terface thickness. At each point r along the interface, 
Eq. pS]) then gives the expression for the interface ten- 
sion 7i2(r) for condensates at effective chemical poten- 
tials /X. — U.{r) (i = 1, 2). The total excess energy ftj^ is 
then obtained by integration over the interface area A i.e. 

= dr7^2(r). The function ^{£,2/^1, K) can be put 
outside this integral as it is position-independent; indeed, 
both K and ^^/^i uniquely determined by the scat- 
tering lengths and the particle masses (see Eq. 0). This 
remarkable property implies that, first of all, along the 
entire interface in a trap, the interface profiles ■01 and 
are unchanged (determined by Eqs. (jlOp and (jllj) onlv). 
This is an exceptional feature of BEG interfaces [5J| . Sec- 
ondly, it implies that in a trap the local interface tension 



7^2(1") depends on position only through the term P^^ 
as can be seen from Eq. (jl3p . The total excess grand 
potential is then: 



4.^(^2/^1,^) / drP(r)^,(r), 



(14) 



which, using the effective chemical potentials /i^ — (r) 
in Eqs. ([5]) and ([9]), brings us to Eq. ([1]) wherein we fac- 
torize the total excess energy as a product of two distinct 
terms: the first is a function of the microscopic gas pa- 
rameters whereas the last is a position integral of the 
external potential along the interface area A. The calcu- 
lation of the function J-' constitutes the subject of next 
section. 

Expression ([T]) for the total excess energy is our main 
result. Together with a Thomas-Fermi (TP) approxima- 
tion in bulk, it forms the core of a simple method to calcu- 
late the total energy of the trapped system and therefore 
to minimize it. The TP approach neglects the energy con- 
tributions arising from density gradients and thus gives 
rise to sharp interfaces, the locus of which is normally 
taken where two-phase coexistence is satisfied. The lat- 
ter, however, is invalid when the interface is curved as 
appears for example for droplets. Instead, a difference 
in pressures appears along the two-phase boundary as 
expressed by the Laplace equation: 



(15) 



with i?j and i?, the principal radii of curvature of the 
interface. It is clear that the difference in pressures along 
the interface is largest for small droplets. Within a first 
approach, for large and R^, one can take the 7^2 in 
Eq. (fT5|) to be the interface tension for a flat interface, 
calculated by assuming equal pressures on opposite sides 
of the interface. 

Pinally, note that, pertaining to the interface tension 
in a trapped system, it is essential to work at fixed chem- 
ical potentials instead of particle numbers and use def- 
inition p3)) as opposed to the interface tension in the 
canonical ensemble [23, ^91 . 



IV. CALCULATION OF INTERFACE TENSION 

In the following four subsections A-D, we derive ex- 
pressions for ^{£,2/^1, K) = 7i2/(4P^i) as expansions for 
different regimes of ^2/Ci and K. The wave functions V'l 
and -02 niust satisfy the boundary conditions (|lip . The 
regions of validity of the expansions are outlined (qual- 
itatively) in Pig. [1] where A-D indicate the subsections 
and the shaded region indicates the absence of an accu- 
rate approximation. In Pig. [21 we depict typical profiles 
for the wave functions "0^ and which are met in the 
four considered regimes. Pirst of all, in A, we focus on 
the case of weak segregation i.e. K ^ 1. Secondly in 
B, we develop an expansion around the point of strong 
segregation 1/K ^ 0; the result is expansion ^ which is 
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1/K 




FIG. 1: Regions of validity of the approximations presented in 
Sects. HVl A-D in the space of parameters 1/K and ^j/Ci- In 
the shaded region, the approximations have a low accuracy; 
the boundaries are drawn on a qualitative basis. The crosses 
denote values corresponding to the experiments of Ref. [13] . 



accurate for a broad range of values. Thirdly, we treat the 
case of a strong healing length asymmetry, that is, when 
^ 1, in C. Finally, we study the case of ^ 1 
in combination with strong segregation 1/ K ^ in D. 



A. Weak Segregation 

For the regime of weak segregation i.e. close to the 
point where the two phases tend to mix [K 1), it was 
found by Barankov in Ref. that the interface tension 
varies as: 



(16) 



In case of ^^/^i — I1 the same leading behavior was re- 
covered in Refs. [H, IH and [H and the square-root be- 
havior as a function of — 1 arises from simple models 
as given in Refs. i21, and [13. Approximation (fT6|) is ex- 
act at K — 1 but is already 3.6% off at if = 1.1 when 
^2 = . Despite this small range of validity, such small 
values for K —1 were realized in the earlyexperiments on 
phase-segregated ®^Rb mixtures [1,[EHI3- the other 
hand, the use of an interface tension for determining the 
excess energy in these experiments is inaccurate as the 
"interface thickness", that is (^^j — 1, is of the same 
order as the length of variation of the trapping poten- 
tial; this invalidates the assumption of constant chemical 
potentials across the interface. 



B. Strong Segregation 

In the following we develop an expansion around the 
point of strong segregation where 1/if 0. An analo- 
gous expansion which originated in a paper by Ginzburg 
and Landau [i^, was used to determine an accurate an- 
alytical expression for the interface tension of a normal- 
superconducting interface [3]. 

In the limit XjK 0, species 1 and 2 do not over- 
lap. Their wave functions are easily found using the 
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FIG. 2: Typical interface profiles for the four cases which 
are treated in Sects. IIVR -D. We indicate the characteristic 
length scales over which the wave functions vary. A) Weak 
segregation regime K ^ \ (here with ^j/Ci = !)• The rel- 
evant length scale is ^^j^/K — 1 which is large. B) Strong 
segregation regime with IjK — > (here with K = 200 and 
^2/^1 = 1). The condensates overlap over a length scale 
(see Eq. (fTT)) ). C) The regime of strong healing length asym- 
metry ^j/Ci 0- When ^j/Ci is small but nonzero, il>r^ van- 
ishes with a tail of length A (see Eq. pS]) ). D) When both 
^2/^1 ~* and XjK — > 0, the condensates will overlap over a 
distance ^2- 



GP Eqs. pn]) and the boundary conditions (|lip : = 
e(2)tanh[z/V2CJ and V2 = e(-z) tanh[-z/V2C2], 
where Q is the Heaviside function. For a finite but small 
value oil/K, these tanh profiles shift towards each other 
so as to overlap over a length proportional to ^/^^/ ^Tk. 
In this overlap zone, the wave functions are modified due 
to the interspecies interactions. In Fig. [2J3, we depict 
such wave functions for the case of ^2/^1 = 1- Note that 
the length of overlap vanishes very slowly when 1 jK — > 
as it is proportional to 1/ \fK. We introduce the follow- 
ing (change of) variables {i — 1,2): 



'K 



and = q„ / 



(17) 



Note that has the dimension of length (drawn in 
Fig. [2J3) whereas and ^2 are dimensionless. We ex- 
pand the wave function of phase i = 1, 2 in terms of the 
parameter <j. : 



We remark that {ip^ — V'")/^; is a regular function of c:'^. 
The functions ip^^ and if)'-^^ are derivable from the asymp- 
totic behavior of ip^ for ^ 00 which are shifted tanh 
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profiles: 



X tanh 



^/2 



^.r,, + ^^r,, + ■ ■ ■ 



X tanh 



V2 



(18a) 



(18b) 



From this expansion, it is clear that ■0°^^ and are first- 
order and third-order polynomials in z respectively. Note 
that the "absolute shift" can be set to zero at every order 
without any loss of generality; thus S-^^ = S^^ and ^'^S^^ = 

Substitution of the expanded wave functions in the GP 
Eqs. (jlOp entails four new differential equations: 



"022 = -^^21 +V'22V'fi + 2-02lV'llV'l2[4/Cl]^- 



(19a) 
(19b) 
(19c) 
(19d) 



Here the dot denotes the derivative with respect to z. 
The solutions of these equations will provide us with the 
numerical values for the (5's; indeed, the boundary condi- 
tions for ip.-^ and 4'i2 to tangentially approach ■0°^^ and 
■0°^ respectively for z — > ±oo, and to vanish for z —^ =Foo. 
Now, in the same fashion as the wave functions, one can 
expand the interface tension: 



7l2 



i=l, 2 



The zeroth-order term is obtained using (unshifted) tanh 
profiles such that 7" = 4V2F (Ci + C2 ) /3 [IS HI , while: 



4^P 



12 1 



0.5 





= 1/2 


^^^^ • ^ 


=1/3 




K=1.5 






— Eq.(3) 




— - Eq.(16) 






M 


0.5 


1 



-1/4 



1 K 



FIG. 3: Values of T or interface tension in units of 4P^j^, 
against values of K~^^'^ where K = G^^j 'JG^^G^^ . The 
squares, triangles and dots are obtained by numerical inte- 
gration of the GP equations for ^j/Ci = li 1/2 and 1/3 re- 
spectively. The full lines are plots of Eq. (|3| and are seen to 
give an accurate approximation for values of K in the interval 
[1.5, oo[. The grey lines plot Eq. ([16} which is accurate for 
T^G [1,1.1]. 



numerical values for the (5's from the numerical integra- 
tion of Eqs. HH). We find that (5^^ = = 0.771 141, 
and, assuming = 1 in Eqs. (|19cP and (|19dp . we ob- 

tain (5^2 — (522 — 0.055. Note that this result is different 
from (^2 found in Ref. H because the problem at hand 
is different starting from the first order in By fitting 
numerically-obtained values for 7^2 with an additional 
fitting term of fifth order in Xj \fK^ and neglecting the 
dependence of and on /Ci j we finally arrive at 
expression ([3]). 

The accuracy of this expansion is clear from Fig. [3] 
where we plot Eq. ([3]) in full lines against the numerically 
obtained values for ^2/^1 = 1 (squares), ^3/^1 = 1/2 
(triangles) and ^2/Ci = 1/3 (dots). Our expansion works 
very well for K in the interval [1.5, +oo[; moreover, it 
is good for all three values of ■^2/^11 which a posteriori 
justifies the neglect of the dependence of , b^^ and the 
fitted higher-order term on ^2 /Ci • Also in Fig. [31 we draw 
with grey lines Barankov's approximation (jl6p . accurate 
for values of K between 1 and 1.1. 



E2 



[4/fi]'0ii'0i2 + 



V'21'022 



After long calculations, one can express the interface ten- 
sion in terms of the spatial shifts (5^^, ^ , ^ , J .,, and ^321 in 
a way similar to what was done in Ref. |44| : 



7„ = 4^, 



3 m 

^^UiAK2\i2li.\ + K2\iJL^) 



It is readily seen in the expansion that 7^2 is symmetrical 
in and ^ it must be. As a final step, we extracted 



C. Strong Healing Length Asymmetry 

So far, we have constructed expansions around extreme 
values for the interspecies interaction parameter K . In 
the following, we will focus on the case when the healing 
length of species 2 is much smaller than the one of species 
1 i.e. when ^2/Ci ^ 1- 

In the limit of ^^2/^1 = (and such that ^2^2 ~^ 0), 
the GP Eqs. (fTO| can be reduced to 

(20a) 



02=e(-z) [i-if^f], 
^f0, =e(z)0, [-1 + 0?] 

+ e(-z)v,(i^-i) [i-(i 



(20b) 
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from which it follows that the density of phase 2 is slaved 
to the density of phase 1. The solution to Eq. (j20b[) 
is (H: 



= e(z)tanh ( 



(21) 



Q{-z)^/2 

VkTT 



cosh 



The constants 2;^ > and < are determined by the 
continuity of "0^ at 2; = where "0^(0) — 1/K. The wave 
function is nonzero only when z < and relates to 
by Eq. (|20a|) . Near z = it vanishes with the square-root 
behavior ip2 oc ■ 

The first-order corrections due to a small and nonzero 
^2/^1 affect only ip^, near the origin z — where is 
large due to the square-root behavior of -02 ■ These cor- 
rections will induce wave function ip2 to acquire a small 
tail of length A which is depicted Fig. We search 
now how GP Eq. (jlObp can be modified to appropriately 
describe such corrections. We assume that "01 varies lit- 
tle on the length scale A; this allows to expand -ft'0j 
in Eq. (jlObp around z = as Kip"^ « 1 + Fz/^^ where 
F = 2 1 K{K — 1). Sufficiently far from the origin, when 

A <C — z <C ^j^, we want to vary as -02 oc z/^^; 
therefore, we may not throw away the nonlinear term -0^ 
in Eq. (jlObp . By appropriately scaling ■02 and z, we find 
that, for |z| < 



= V^2(^ + V'^)' 



where the dot denotes the derivative with respect to z 
which we define as: 



z/A with \ = iAF\Ui.\\ 



-1/3 



(23) 



(22) 



and where the new wave function i/'^ = "02 [^[^2/^1]] ^^'^ ■ 
For — z ^ A, the term -0^ in Eq. may be neglected; 
indeed, doing so results in the solution -0^ — \/—'z. As 
expected this is in fact ijj^ oc — y^/^i- 

The differential equation which appears in Eq. ((22|) 
was encountered earlier by Lundh et al. [32] and Dalfovo 
et al. [m when studying the BEG order parameter at 
the border of a large harmonic trap. The role of the 
trapping potential is played here by the interactions with 
condensate 1. 

We briefiy sketch how we calculated the interface ten- 
sion; we refer to Refs.[3lland[3^! for details of the method. 
To zeroth order in /^i i the interface tension is easily 
found by integration of -0^ (see Eq. (1131) ) with 0^ given in 
Eq. (j2ip . For small and nonzero ^2 /^i > extra energy con- 
tributions arise due to the integral over 02 (see Eq. ([13])). 
For — z ^ A, -02 must be taken from Eqs. (|20ap and ((2T|) 
whereas for |z| ^ 02 is the (numerical) solution of 
Eq. (|22p. Since these two solutions for -ij)^ are equal in 
the region A ^ — z ^ , where V'a oc a/— z/^j the inter- 
face tension is well-defined. All the extra energy contri- 
butions due to a nonzero ^2 /^i can be seen to be of order 

K2/CJ'. 

Galculations then lead to T to second order in ^2 /^i • 



7l2 




K2/e. 



,2 V^£T ( K + i 



— arctanh 



IK - 1 
2K 



(24) 



Since there are two length scales (A and ) over which the 
density varies, their relation A <C renders difficult the 
full numerical integration of the GP Eqs. Accordingly, 
we could not verify expansion (|24p directly. 

Gonditions for the validity of the above approxima- 
tion can be derived. First, the linearization of 0^ near 
z = in Eq. (jlObp is justified only when the length 
over which higher-order corrections to 0"^ are relevant 
(that is [i/'^/V'J(0) oc V^^J is large as compared to the 
length over which ip^ varies (being A). Second, correc- 
tions to 01 caused by a nonzero ^2/Ci can be neglected 
on [—A, 00 [ when, in Eq. pOap . the amplitude ^^"02 of the 
term Kijj^il)^ (being ii'[F[^2/^i]]^^'^) is sufficiently small. 



These conditions of validity amount to: 



« ^K{K-l) and ^2/^1 « Vif. 



Thus, our expansion fails in the regime of strong segre- 
gation when ^/K ^ £,2/^1 and close to weak segregation 
when ^JK — \ <C ^2/^1- Note that to zeroth order in 
5,2/^11 our expansion (|24|) agrees with (fTG]) when K 1. 
However, it differs from the expansion (fTB)) by a factor of 
2 in the term of order [Ca/CJ^- 
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D. Both Strong Segregation and Strong Healing 
Length Asymmetry 

One may now wonder what happens to the interface 
tension when both K and [^2/?!]""'^ ^'^s large. When K 
goes faster to infinity than [^2/^1]"^: '^^ will find that 
expansion ([3]) is still valid while expansion (|24p must 
be taken in the inverse case. In the following, we fo- 
cus on the intermediate case when both [■^2/^1]^^ ^^id 
K diverge such that k = ([^2/Ci]v^)~^ remains of order 
one. To zeroth order in both ^2/^1 the conden- 

sate wave functions do not overlap; seen on a lengthscale 

, this results in the wave functions = Q{—z) and 
V-i = e(z)tanh[z/V2Ci]. 

For small and nonzero values for / Ci 1/K while 
AC = {[£,2/(.i]VK)~^ is of order one, the condensates over- 
lap. Due to the strong repulsion this happens only within 
a short region of length ^, ; over that interval, the value of 
the wave function varies between zero and one while 
the value of ip^ does not vary much (see Fig. [2P). The 
latter is only modified close to z = 0, where, to zeroth or- 
der, it vanishes linearly. This brings about the definition 
of a new wave function 0, : 



result, one finds: 



4P^, 3 ^^^/^^M 3 ([4/^jVx)i/2 



7l2 



(27) 



0.06653 



0.00107 



This approximation is very accurate for values of k in the 
interval between and 1 44]. To first order in k, Eq. ([27|) 
agrees with the interface tension found for K ^ 00 as 
given in Eq. ([3]). 

For large k, according to standard results, the integral 
in Eq. ^ goes as 2(1 - \/2)k/3 such that 



7l2 

4PCi 



V2 , 2(1 - ^/2) 



3VK 



+ [^J^^]0{n-'). (28) 



Since large k implies ^ vK, this last result is in 

full accord with the K 00 limit in expansion Eq. (j24p . 
The higher-order corrections to Eq. (|28p can also be ex- 
tracted from Eq. 

Note that Eqs. (pS)) and are valid for all values of 
K and may be generally used to numerically integrate the 
interface tension in the limit under consideration. 



V', = ^2/^,] 



z'e(z') 



V2 



Q(z) tanh 



(V2C, 



Here we introduced the dimensionless variable z' = z/^^ 
and (j)^ must have the asymptotic behavior 0^ (z' ^ 00) = 
z' / a/2 and 0^ (z' —00) = 0. We can now expand the 
GP Eqs. pn]) in terms of the small parameters [^2/^1]^ 
and 1/K in the overlap interval: 



(Pj_ \ K / <f'2 \ '"^ 



(25) 



The dots denote the derivative with respect to z'. One 
may notice that these equations are the same as the 
Ginzburg-Landau equations which are valid in a super- 
conductor when a constant magnetic field is applied par- 
allel to its surface. In that case, the Ginzburg-Landau 
parameter k is the ratio of the penetration length of the 
magnetic field to the coherence length, tp^ plays the role 
of the order parameter of the superconductor and 0^ that 
of the vector potential [5^1 . 

By a straightforward calculation, one can rewrite the 
interface tension Eq. (|13p as an integral over 0^ and : 



7l2 

4e,p 



dz' 



1 



(26) 



where the dots denote the derivative with respect to 
z' and and ip^ obtained by the Euler-Lagrange 
Eqs. (ESI). 

An analytical expression for the interface tension of a 
superconductor-normal interface in case of small k was 
obtained by Boulter and Indekeu in Ref. 44. Using their 



V. TEMPERATURE DEPENDENCE OF THE 
INTERFACE TENSION 

At low but nonzero temperature, interface waves will 
be excited as thermal excitations. In the following sec- 
tion, we will consider these and ignore quantum fluctua- 
tions. It is common to incorporate the free energy contri- 
butions which arise from capillary excitations into a new 
temperature-dependent interface tension j-^^{T). Gener- 
ally, when the interface waves have wavenumbers k and 
frequencies uj(k), one finds for a three-dimensional sys- 
tem that f47ll48ll: 



7i2(r) = 712(0) + 



27r 



In 1 



fcdfc. 



where 7^2 (0) is the interface tension at zero temperature 
and the last term embodies a negative correction due 
to a nonzero temperature. For a flat interface between 
two BECs, one may obtain the capillary-wave spectrum 
from the time-dependent GP Eqs. For sufficiently large 
wavelengths ^ 712/^1 '^^^ then finds for the modes 
with in-phase motion of both species [1^ [13] : 



7l2 



(29) 



where rii is the bulk density of phase i — 1, 2. 
When the temperature T satisfies the condition k^T <^ 
fi^/(2m.(7j2/-P)^), only modes with wavenumbers satis- 
fying Eq. (|29p are thermally excited. The temperature- 
dependent interface tension then becomes ^56^]: 



7i2(r)- 712(0)- 



47rfi4/3 



2/3 



712(0) 
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C3 

85DK 



U-{r,z) = m^{uj1r'^ + a;^.z^)/2 with = a;^ + y^. Us- 
ing Eq. ([T|), one finds the total excess energy due to the 
presence of the two interfaces: 



FIG. 4: Cross section of a trapped cloud of BEC mixtures as it 
was observed in Ref.fiol. Due to phase segregation, *"'Rb forms 
a layer between clouds of ^'^Rb. In Eq. (|30|l . we calculate the 
total excess energy due to the presence of the two interfaces. 



where c = r(7/2)C(7/2) and C, and T are the zeta and 
gamma function, respectively. 

Taking ~ n^, 712(0) oc and /i^ ~ /^j, one calcu- 
lates that the relative temperature corrections [j-^^ (T) — 
iMh.M are of order ^^^^^(fc^T//. JV3_ 

known that ^n^a\^ ^ 1 is the condition for the GP 
theory to be applicable; moreover for the current BEC 
experiments T is of the same order or less than [4ll | . 

On the other hand, for a system in the weakly seg- 
regated regime, the interface tension varies as (0) cx 
£,^P^/K^^ and accordingly vanishes as i^T — > 1. The 
relative temperature corrections are then of order {K — 
\)~^/^ ^JrLjjt^^{kgT / /'^ which may become large as 
K ^l. 

In conclusion, thermal fluctuations of a flat BEC in- 
terface do not affect the interface tension at low temper- 
ature, except possibly near weak segregation. Note that 
the calculated temperature dependence of the interface 
tension does not hold in trapped systems, as the inter- 
face modes there are quantized [13, HH. 



VI. DISCUSSION 

We discuss now the relevance of our results with re- 
spect to the recent experiments of Papp et al. p^ . 
Two main values for the scattering lengths of phase- 
segregated states were reported: one configuration was 
characterized hy K = 3.01 and ^2/"^! — 0.84 such that 
T — 7i2/(4P^J = 0.434, while for the other configura- 
tion, K = 2.36 and = 0.952 such that T = 0.415. 
Here species 1 and 2 are ®^Rb and *^Rb respectively. 
Note that these values for T are obtained both numeri- 
cally and with Eq. ([3]). The values for T together with 
Eq. ([T]) allow a numerical calculation of the interface ex- 
cess of any trapped configuration. 

As an example, we consider now the topology as de- 
picted in Fig. [3] in which a layer of ^^Rb BEC is present 
between clouds of *^Rb. This configuration was experi- 
mentally obtained in Ref. [l^ by strong radial confinement 
so as to quench gravitational effects. Particles of phase 
i = 1, 2 are confined by an anisotropic trapping potential 



QHO 



m^a;;^^(r° -I- r°) 
20fia,, ' 



(30) 



Here and are the maximal radial coordinates along 
the two interfaces (see Fig. 0]) and T — 0.415. Com- 
paring this excess energy with the total energy 
i? of a gas of 1.4 x 10^ ^^Rb particles, we obtain that 
/E Ri 0.03. This irnplies that, even for the largest 
clouds observed in Ref. [lO, the interface excess energy 
is substantial and must be considered accurately when 
determining the shape of the ground state configuration. 

Note that in recent experiments on ultracold imbal- 
anced fermion gases, a breaking of the local density ap- 
proximation was encountered in the sense that the intro- 
duction of the interface and its tension was essential in 
explaining the experiments [iol [sot . 

The formalism presented in this work assumes a lo- 
cal approximation for the interface tension. For this 
approach to be valid, the interface thickness must be 
smaller than the length over which the trapping po- 
tential varies. Generally, the former length is simply 

-1-^2! however, close to the mixed state (i.e. when 
K ~i 1) it is (Ci -h i^)l\/K - 1 which is larger. On 
the other hand, the characteristic variation length o f the 
trapping potential of species i is a^^ ~ y^h/m.tJ. with 
uJ = {lu^ .Lu^ .y^^ . If we consider the configuration of 
1.4 X 10^ particles with K = 2.36 and ^2/^1 = 0.952, 
we find that (^^ + £.2)/(^ho ~ 0.1 is indeed small at the 
trap center [lOj . Although the coherence lengths diverge 
at the trap boundary, this divergence is so weak that 
Ci + ^2 exceeds only in a layer of thickness of 0.5% 
of the trapping radius. This justifies the use of the lo- 
cal approximation for the interface in the experiments of 
Ref. m 



VII. CONCLUSION 

We propose a method to calculate the total energy of a 
trapped phase-segregated BEC mixture. In order to find 
the bulk energies, one can use a local density approx- 
imation or Thomas-Fermi approach. However, with re- 
spect to the recent experiments on phase-segregated BEC 
components [l^], the bulk contribution must be supple- 
mented by interface excesses which may have pronounced 
effects on the ground state topology. The interface excess 
energy arises from the presence of the two-phase bound- 
aries, the position of which are determined by the Laplace 
equation. Using a local approximation for the interface, 
we write the total excess energy as a simple integral of the 
trapping potential along the interface (see Eq. ([T|)). By 
use of analytical expansions, we find the interface tension 
in almost all parameter regimes (see Fig.[T]), including the 
regimes encountered in Ref. [l^. We also study the influ- 
ence of thermal fluctuations on the interface tension at 
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low temperature and find that their effect is negligible, 
except possibly close to the mixed state. 
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